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The Banks-Casher relation links the spectral density of the Dirac operator with the existence of a 
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lattice, from the finite value of the lattice spacing a. We present a status report of a determination 
of these corrections for Wilson quarks. 
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1. Introduction 

The Banks-Casher relation [1] relates the spectral density p/j of the Hermitian Dirac operator —iD 
with the chiral condensate £ 

PD{y,m) = ^[l + 0{{\y\,m]/^QCD)] , (1.1) 

where are the eigenvalues of the massless operator and m is the sea quark mass. In principle the 
Banks-Casher relation provides a tool to determine the chiral condensate [2] with a lattice QCD 
computation. 

A more standai^d way to determine the chiral condensate with lattice QCD computations is 
to study the quark mass dependence of the pion mass and comparing it with the predictions of 
chiral perturbation theory (;tPT). Recent lattice calculations are performed close to the physical 
values of the quark masses (see ref. [3, 4] for recent reviews). To control accurately the light quark 
mass dependence of hadronic quantities it is important, if possible, to have independent determi- 
nations of leading order (LO) low energy constants (LECs), as for example the chiral condensate. 
In the plot in fig. 1 we show the light quark 
mass dependence of the squared pion mass com- 
puted by the European Twisted Mass Collab- 
oration (ETMC) with Nf = 2 dynamical light 
quarks [5]. Despite the rather good description 
of the lattice results with ;^PT, it is obviously de- 
sirable being able to predict, in an independent 
way, the behaviour close to the chiral limit (red 
ellypse) to better constrain the light quark mass 
dependence. A constraint of the chiral fit will be 
beneficial for a more accurate determination of 
the next-to-leading order LECs and a better con- 
fidence on the chiral fits. 

An example of independent determination of 
LO LECs is provided by lattice calculations in 
the so called epsilon-regime, where first estimates 
of the chiral condensate and decay constant are 
rather encouraging [6, 7]. Using Wilson-type fermions as a QCD discretization, calculations are 
affected by potentially large cutoff effects, thus it is important to have a theoretical analysis of quan- 
tities such as the spectral density, based on Wilson chiral perturbation theory (W;^PT) [10, 11]. ^ 
These proceedings are a status report of an ongoing attempt to understand cutoff and finite size 
effects affecting the spectral density of the Wilson operator. 

2. Chiral condensate from the mode number 

With Wilson fermions it is advantageous to consider the Hermitean Wilson-Dirac operator Q = 
y^Dm, where D,„ is the massive Wilson opearator. To compute the spectral density pg in W;^PT 

'The analysis of two-point functions in the epsilon-regime within the framework of Wj^jPT has been discussed in 
refs. [12, 13] and more recently extended in ref. [14]. 
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Figure 1: Light quark mass dependence of 
the ratio between the pseudoscalar meson mass 
squared and the renormalized light quark mass 
in units of the Sommer parameter tq. 
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Figure 2: Left plot: the spectral density [pg(A,m) +Pg(— A,m)]iVLo in the infinite volume. We used the 
parameters I = (275 MeV)^ otpcac = 26.5 MeV, F = 90 MeV, Lg = 5, = 139.6 MeV. The solid black 
line corresponds to the continuum ;i;PT prediction, while the red (blue) lines correspond to the lattice W;);PT 
prediction (6'(a)-improved) on eq. 2.3 with a^Wg = ±5 • 10^ MeV*, corresponding to A = =f0.53. Right 
plot: result of the global fit of the data published in [2] with our formula for the renormalised mode number 
(eqs. 2.4 and 2.3). The fit parameters are E, A and Lg and we obtain ;|;^/dof = 0.91. 

one introduces a flavour doublet of valence fermions Xv with a Wilson twisted mass action and 
twisted mass ji^,. The spectral density is related to the discontinuity of the valence pseudoscalar 
condensate along the imaginary axis in the twised mass plane [15] 



Disc [{XvYs^^Xv, 



lin [pq{X ,m) + Pq{-X ,m)] 



(2.1) 



thus the valence pseudoscalar condensate is a tool to compute the spectral density in W;^PT ^. After 
matching the continuum Symanzik effective theory with the generalized effective chiral Lagrangian 
one needs to choose a proper power counting for the scales involved in the problem. The scales are 
the sea quark mass, m, the valence twisted mass pLy (directly related to A), the lattice spacing a and 
the linear size of the space-time volume L. In the following we always consider the sea quark mass 
in the so-called p-regime. For the other scales we consider m ~ /Xv ~ a ~ l/Lr^O{p^). Theresuh 
of the calculation [16] in terms of the PC AC quark mass mpcAC is given by 



[pQ(A,mpcAc) +Pe(-'^,'"PCAc)]wLO = 2[pQ(A,nipcAc)] 
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(2.2) 
(2.3) 
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^ + ^ + ^ ) > Mvs is the pseudoscalar meson mass (made of two sea 



In these formula A 

quarks), £ and F are the LO LECs, a = IWoa and Wo as the others Ws are the LECs parametrizing 
0(a) and O(a^) effects. Details on the calculation and a discussion on the applicability of this 
formula can be found in refs. [16, 17]. 

Potentially even with Wilson fermions one can use the spectral density, or equivalently the 
renormalization group invariant (RGI) mode number [2] 



.A 

v{A,m) = V dXpQ{X,m 

J-A 



(2.4) 



^We recall that in the continuum the two spectral densities are connected by the relation pQ[X,m) = 
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to determine the chiral condensate. To test our formula we compared the resuhing mode num- 
ber (2.4) with the numerical data published in [2]. We have fixed F = 90 MeV and the renormali- 
sation scale }X = m^; we have performed a global fit at all the 3 masses available and all the values 
of Ar with 3 fit parameters: £, A and L^.^ From the global fit we obtain 

1^/^ = 266(7) MeV, A = -0.62(80) , L6 = 6(l). (2.5) 

The numerical data and our global fit are shown in the right plot of fig. 2. For £ we obtain a 
perfectly consistent result with Giusti and Liischer [2] without performing any chiral extrapolation. 
We have also performed a fit fixing F = 80 MeV obtaining compatible results within errors. Lattice 
determinations of the chiral condensate using the mode number can be found in [2, 8, 9]. 

3. Cutoff effects close to the threshold 

To improve the theoretical description of the spectral density of the Wilson operator for A ~ m 
one needs to consider two important points. The first one is that finite size effects diverge for 
X ^ m [16] and the second one is that the when A ~ m the power counting for the valence quark 
masses need to be reconsidered, i.e. the treatment of the cutoff effects in a perturbative fashion 
might not be adequate. 

To overcome this difficulties we opt for the following power counting 



m^O{p'), mp=^m2+ju2~0(p"), l/L,l/T ^ 0{p) , a ^ 0{p') , (3.1) 

that implies cutoff effects affecting the spectral density at NLO order. The framework is the so- 
called mixed Chiral Effective Theory [18], where some masses obey the /7-regime counting, and 
others are in the epsilon-regime. As an intermediate step of our calculation we introduce a 0-term 
solely in the mass-term of the action as follows 



F 
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^2 = ^Tr [d^Ud^U^ - -Tr \ulu{x)^^ + ^'^U{x)U0'^ - — Tr [U + U^ , (3.2) 
where the mass matrix in the replica formalism [19] is 

Ji = Jt'^ = diag(m, ...,m,mv + //ivT^,---,'Wv + 'MvT^), f/e = diag(e^"'^% Ir) . (3.3) 

While in the continuum it is not important how the 0-term is introduced in the parametrization 
of the {/-field [18, 20], adding a 0-term in the sea sector only, becomes relevant at finite lattice 
spacing. With this choice, even with a 0-term in the action, we can reabsorb the leading O(fl') 
cutoff effects in a redifinition of the quark mass, because with our choice of power counting (3.1), 
in the sea sector the leading O(fl') effects appear at NNLO, i.e. the sea quarks are effectively in 
the continuum up to higher order corrections. With this particular- choice of power counting and 
parametrization of the 0-term we can achieve, as in the continuum [18, 20], a factorization of the 
partition function for the zero and non-zero modes. The periodicity in Q of the chiral Lagrangian 
allows us to write the partition function in standard fashion 

3^(6)= y e-"'^^, irv = — / dde'^^^ie). (3.4) 
vi^oc 271 Jo 



^In this proceedings with Z we denote its value renormalised in the MS scheme at a scale of 2 GeV. 
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By performing an exact integration over the constant field Q one obtains 



dUo (detUo) 



Tr U'o+Uq 



Zs = mVL , 



(3.5) 



from wliicli one observes that the distribution of v is Gaussian and it is controlled by the sea quarks 
which are in the p-regime. The computation of the spectral density is now a computation at fixed v, 
i.e. p = £v Pv The sum over v can be done because we know the weight factor (cfr. eq. 3.5). 
We expand the pseudo Nambu-Goldstone field U (x) around the ground state of the theory Uy 



J7v=diag( 1 ,e sma)o = — , coscoo = — 

N, N,. 



(3.6) 



and with this parametrization, the mass term in the chiral Lagrangian becomes like in the untwisted 
case, with a degenerate polar mass nip in the valence (replicated) sector. 

We compute the pseudoscalar valence condensate (cf. eq. (2.1)) in the chiral effective theory, 



where now the source term has the following form 



The final result of the calculation can be written as 

^eff , 



17=0, 



(3.7) 



1^ 



:diag( , 0,...,0). (3.8) 



Nr 



+ - [ABdW(,NsrriVF2 + d^VW^{F^ + AN,F2) + d^VWjFs+a^VW^F(] , 



(3.9) 



where 7^=i..6 are integrals over zero modes and Zeff is an effective chiral condensate 
which includes NLO corrections from the 
sea quarks and the lattice spacing 



^eff=(reff)cont + 16S^. (3.10) 

The explicit formula for (reff)cont ^^"^ 
found in [21]. In fig. (3) we show the impact 
of the NLO corrections to the chiral conden- 
sate [21]. One observes that in the typical 
range of quark masses for dynamical simula- 
tions the relative corrections induced by the 
presence of sea quarks can reach 30%. This 
is a warning in case one would like to extract 
the chiral condensate from fits of mode num- 
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Figure 3: Significance of the NLO corrections for 
Eeff/E in a typical range of sea quark masses for dy- 
namical lattice QCD computations. 



bers where the effect of the sea quarks have been neglected. To compute the zero-modes integrals 
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we switch to the supersymmetric formulation. All integrals at fixed v can be computed by derivat- 
ing the graded SU {2\2) partition function [22, 23] with respect to appropriately chosen sources. 
The integrals have been computed and cross-checked. The final step to determine the formula for 
the spectral density is the calculation of the discontinuity along the imaginary axis in the twisted 
valence mass plane (cfr. (2.1)). The analysis of the final result is in progress. We conclude ob- 
serving that the effects of the sea quarks are twofold. They change the absolute normalization by 
introducing an m dependence in Zgff and they control the distribution of v, which, we remark, stays 
Gaussian only because the sea quarks are in the /^-regime. 
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